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An antiferro ordering of an electric hexadecapole moment is discussed as a promising candidate for the 
long standing mystery of the hidden order phase in URu 2 Si2. Based on localized /-electron picture, we 
discuss the rationale of the selected multipole and the consequences of the antiferro hexadecapole order of 
xy(x 2 - y 2 ) symmetry. The mean-field solutions and the collective excitations from them explain reason- 
ably significant experimental observations: the strong anisotropy in the magnetic susceptibility, character- 
istic behavior of pressure versus magnetic field or temperature phase diagrams, disappearance of inelastic 
neutron-scattering intensity out of the hidden order phase, and insensitiveness of the NQR frequency at 
Ru-sites upon ordering. A consistency with the strong anisotropy in the magnetic responses excludes all 
the multipoles in two-dimensional representations, such as (O yz , O a ), The expected azimuthal angle depen- 
dences of the resonant X-ray scattering amplitude are given. The (x 2 -y 2 )-type antiferro quadrupole should 
be induced by an in-plane magnetic field along [110], which is reflected in the thermal expansion and the 
elastic constant of the transverse (cu - Cn)/2 mode. The (x 2 - y 2 )-type [(xy)-type] antiferro quadrupole is 
also induced by applying the uniaxial stress along [110] direction [[100] direction]. A detection of these 
induced antiferro quadrupoles under the in-plane magnetic field or the uniaxial stress using the resonant 
X-ray scattering provides a direct redundant test for the proposed order parameter. 

KEYWORDS: multipole order, neutron scattering, resonant X-ray scattering, Holstein-Primakoff method 
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1. Introduction 

The heavy-fermion superconductor URu^Sia 1-3 ' undergoes 
a phase transition into the mysterious phase referred to as 
"Hidden Order" (HO) at T a = 17.5 K at ambient pressure. The 
term, HO, is given by the fact that the phase transition is sig- 
naled clearly by a large anomaly in the specific heat, while its 
order parameter (OP) has not been identified by a number of 
experimental investigations over the last two decades. In early 
stage of research, the observed tiny antiferro (AF) magnetic 
moment 4 ' 5 ' was considered as a key factor for the HO phase. 
However, it turned out that it is the parasitic tiny moment, 
whose volume fraction is enhanced by applied pressure. 6 ' 7 ' 
The elastic neutron-scattering experiments underpressure re- 
vealed a sharp phase transition at p c ~ 1.5 GPa, above which a 
type-I antiferromagnetic (AFM) phase emerges with 10-times 
increase of the magnetic moments (/i ~ 0.4//b) along c axis 
(the ordering vector is Q — (100)). 8 ' In contrast to the dras- 
tic increase of the magnetic moment across the /^-induced 
phase transition, the /^-dependence of the transition temper- 
atures from the paramagnetic phase to the low-/? HO and the 
high-/? AFM phases is weak and continuous 9-1 1 > (the transi- 
tion between two ordered phases is of the first order). Note 
that the bulk superconductivity below T c ~ 1 .5 K does not co- 
exist with the high-/? AFM, 9 ' although similarity of the Fermi 
surface in both the HO and the AFM phases is pointed out by 
Shubnikov-de Haas measurements. 12 ' The recent laser angle- 
resolved photoemission spectroscopy gives an evidence of the 
broken translational invariance in the HO phase. 13 ' 

Meanwhile, it has developed that multipole moments of the 
localized / electrons are regarded as a key concept to under- 
stand HO phases observed in various /-electron systems. 14 ' 
This concept should be useful when a long-range order is trig- 
gered predominantly by the localized component of / elec- 
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trons, even though URu2Si2 is one of the most spectacular 
examples showing the itinerant-localized duality inherent in 
the heavy fermions. 15-18 ' So far, the various rank of multi- 
pole orders have been proposed for the HO phase: the mag- 
netic dipole, 19 ' the electric quadrupole, 20 ~ 23 > the magnetic oc- 
tupole, 24 ' the electric hexadecapole 25 ' 26 ' and the magnetic do- 
triacontapole. 27 ' In principle, the resonant X-ray scattering 
(RXS) can detect multipole orders up to the 4th rank, 28 ' 29 ' 
provided that a superlattice reflection from the E2 (electric 
quadrupole) transition has a measurable intensity. The re- 
cent RXS study in a major region of k = (hOl) plane has 
observed no traces of significant intensity except from the 
AFM, indicating the absence of quadrupole orders within 
a searched k range including a conceivable nesting vector 
Q* =(1.4,0, 0). 30) 

Recently, Haule and Kotliar have proposed the hexade- 
capole scenario based on the electronic structure computation 
combined with the dynamical mean-field approximation, 25 ' 26 ' 
where the Kondo renormalization process is arrested by the 
competition with the crystalline electric field (CEF) splitting, 
yielding a self-induced OP across the CEF splitting. Their 
scenario provides a natural explanation for similarity of two 
ordered phases, although they are characterized by distinct 
OPs. It is also pointed out that the low-energy intensity seen 
by the inelastic neutron scattering only in the HO phase can 
be understood by the pseudo-Goldstone mode from almost 
energetically degenerate ordered phases. 26 ' 

Moreover, a space group analysis has been carried out 
to narrow down possible candidates for the OP in the HO 
phase. 32 ' The second order phase transition from I4/mmm 
(No. 139) to P42/mnm (No. 136) does not need any kind 
of lattice distortion, allowing the NQR frequency at Ru-site 
unchanged. On the other hand, the AFH order corresponds 
to No. 128 (P4/mnc), which loses the four-fold axis at Ru- 
site, and it should affect the NQR frequency in contrast to the 
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No. 136 order. However, the asymmetry parameter of the field 
gradient at Ru-site might be quite small due to the higher-rank 
AFH order as will be shown in this paper. Thus, no anomaly 
observed in the Ru-NQR frequency within experimental ac- 
curacy 33 ' gives no decisive choice among two plausible space 
groups. 

Motivated by these circumstances, we first discuss the ra- 
tionale of the selection of multipoles based on the localized 
/-electron picture and the group theoretical point of view. The 
observed strong anisotropy in magnetic responses requires the 
singlet-singlet-singlet low-lying CEF model, among which 
the dipole, the hexadecapole and the dotriacontapole remain 
active. We discuss that the hexadecapole is the most natural 
candidate for the HO, and then we elucidate the consequences 
of the AF hexadecapole (AFH) order of xy(x 2 - y 2 ) symme- 
try in detail. This paper is organized as follows. In §2, we 
classify the CEF basis and the multipoles under tetragonal 
Dm symmetry. We narrow down relevant active multipoles, 
mainly according to consistency with the strong anisotropy 
in the magnetic susceptibility. In §3, we investigate the mean- 
field solutions and the collective excitations from them, which 
explain reasonably significant experimental observations. The 
emphasis is placed on the induced quadrupole moment under 
the in-plane magnetic fields, which could provide a useful test 
for the proposed OR The final section summarizes the paper. 

2. Selection of Hidden Order Parameter Candidate 

In order to narrow down possible candidates of the hid- 
den order parameter, let us classify the multipoles within 
J = 4 multiplet in f 2 configuration. Under the tetragonal 
(Dm) symmetry, which is characterized by the CEF hamil- 
tonian in terms of the Stevens operators, 14 > 34 - 36 ) 

Hqef = B20O20 + B40O40 + B44O44 + BsoOm + BmOm, (1) 

the CEF states are summarized in Table I. Here, instead of the 
five CEF parameters B m „, we use the center-of-mass of ener- 
gies, Ai, A34, A5 and the mixing parameters, 9 (0 < < 7t/2) 
and <p (0 < (f> < n 12) characterizing the CEF wave functions. 
The relation between two sets of parameters is given by 

B20 = -^( 5A i +4A34 + 6A5), 



Table III. The active multipoles between the classified CEF states. 
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where A = 3Ai + 2A 3 4 + 4A 5 and A' = 23Ai + I6A34 + 12A 5 
for notational simplicity. For later convenience, the CEF states 
are classified into three groups, si (+1), 52 (-1) and d (0) by 
the character of C4 operation (as indicated in the parenthesis). 
Note that the groups s\ and d are independent with each other. 

By means of the irreducible representation of the tetragonal 
point group, the multipoles up to 5th rank are defined as in Ta- 
ble II. The corresponding quantum operators are obtained by 
replacing the polynomial of (x, y, z) with the symmetric prod- 
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Table IV. The proposed models and the candidates of the order parameters. 
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uct of (J x , Jy, / z ). 14 ' 34 ' Here, we again classify the multipoles 
into three groups, Si, S2 and D by the character of C4 op- 
eration. This classification is convenient to discuss the active 
multipoles with finite matrix elements among the classified 
CEF states. 

The active multipoles between the classified CEF states are 
summarized in Table III. The multipoles having unit character 
of C4 operation remain active in the "diagonal" parts, si and 
Si, or S2 and 52 in Table III. Similarly, the "off-diagonal" part 
between si and S2 has the active multipoles with -1 charac- 
ter of C4 operation. The multipoles with ±1 character of C4 
operation remain in the diagonal part between d and d. The 
multipole operators belonging to the two-dimensional repre- 
sentation (zero character of C4) appear in the "off-diagonal" 
part between (s\ , S2) and d. 

With these preliminaries, let us examine the so far proposed 
OPs based on the multipole ordering models, which are sum- 
marized in Table IV. According to the analyses of the mag- 
netic susceptibility, x(T), the most natural choice of the low- 
lying CEF scheme is r ( 1 1) -F2, 37) which is common in the re- 
lated materials, UPt 2 Si 2 , 19) PrRu 2 Si 2 38 ' 39) and their dilute sys- 
tems. 40 ' Note that this choice for the low-lying CEF scheme 
is also reached by the electronic structure computation com- 
bined with the dynamical mean-field approximation (Ti and 
T2 are reversed due to the hybridization effect). 25 ' A promi- 
nent feature ofx(T) in URu2Si2 as compared with the related 
materials is the strong anisotropy, namely, Xa(T) shows the 
Van Vleck type almost T-independent behavior, while Xc(T) 
exhibits the Curie- Weiss behavior at least below 300 K. The 
strong anisotropy is observed in the magnetization process as 
well. 19 ' This indicates that the excited Fs doublets are located 
at least 300 K far from the low-lying CEF singlets, which 
makes J x and J y operators inactive at low temperatures. The 
large energy separation between the singlets and the doublets 
results in all the multipoles in the two-dimensional represen- 
tation (the group D) inactive, and hence they can be excluded 
for possible candidates of the hidden order parameter. 

It should be noted that the observed Xc(T) in URu2Si2 
shows a peak structure at about T* ~ 40 K. This is incom- 
patible with the simple analysis with the above Tj 1 '^ CEF 
model, showing the Van Vleck saturation below the F ( 1 1) -F2 
splitting. Practically, the Fermi-liquid quasiparticle emerges 
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Table I. The CEF states of J = 4 multiplet under Dm symmetry. The CEF states are classified by the character of C4 symmetry operation for later 
convenience, which is indicated in the parenthesis of the first column. 



Group Symmetry Wave function Energy 
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Table II. The definition of multipoles up to 5th rank under D4;, symmetry, which are classified by the character of C4 symmetry operation. 
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below T* where the non-trivial inter-site effect should be 
taken into account beyond simple localized models. In fact, 
the decreasing behavior of Xc(T) below T* is diminished by 
Th, La and Y substitutions for U in URu2Si2- 41 ^ 3 ' Moreover, 
the comparable energy scale of T* and the low-lying CEF 
splitting is consistent with the fact that the dilute U alloys 
exhibit non-fermi-liquid behavior, 44 ^ 16 ' which can be under- 
stood by the competition between the Kondo singlet and the 
CEF singlet ground states. 47 ' 49 ' 

The models proposed by Nieuwenhuys, 19 ' and Haule and 
Kotliar 25 ' are consistent with the susceptibility analyses. On 
the other hand, the model by Ohkawa and Shimizu 20 ' could 
explain solely the strong anisotropy in x(T), provided that 
the rest of the CEF states are far from the ground doublet. In 
this case, however, we would have to consider that URu2Si2 
is an exceptional in the series of compounds. Note also that 
no indications of the quadrupole orders predicted by this 
model have been observed especially by using the RXS mea- 
surement. 30 ' 31 ' In the model by Santini and Amoretti, 21-23 ' 
it is assumed the large exchange coupling constant for their 
quadrupole order, O xy or 6*22 and the possible J z order. This 
assumption involves much higher CEF states, which results in 
many active multipoles including J x and J y . Consequently, the 
strong anisotropy inx(T) disappears. The model by Kiss and 
Fazekas 24 ' has a similar difficulty where the low-lying singlets 
and the doublet in their model make J x and J x active. In their 
octupole ordering model, the crossing of the singlet and the 
doublet states is indispensable for their main conclusions in 
the magnetization process and the p-induced transition to the 
AFM phase. 

Eventually, the models by Nieuwenhuys, 19) and Haule and 
Kotliar 25 ' could provide the most plausible candidate for the 
HO in view of the CEF scheme. When we take the group s\ as 
the low-lying CEF states, the multipoles in the group S 1 are 
natural candidates for the hidden OP occurring at low tem- 
peratures. Apart from the possibilities of staggered "scalar" 
(Tig) orders as was observed in PrRu4Pi2, 14 ' the possibilities 
of T2g-type hexadecapole {H") and the Fi„-type dotriacon- 
tapole (D S 4) orders remain (the F2„-type dipole J z is of course 
for the high-p AFM phase). The local wave functions un- 
der J z , H" and D s n orders are schematically shown in Fig. 1, 
where the shape and the color-map represent the charge and 
the magnetic-charge densities, respectively. 34 ' The similarity 
of the local wave functions of Fig. 1(a) and (b) indicates slight 
change of whole electronic states across the ^-induced phase 
boundary between the AFH and the AFM. As was pointed 
out by Haule and Kotliar, 25 ' 26 ' the AFH and the AFM sce- 
nario explains other prominent experimental observations. 26 ' 
In the next section, we discuss the consequences of these two 
orders in detail. 

3. AF Hexadecapole Order and its Consequences 
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The phase diagram based on T^P-Tx-T^ CEF model 
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Let us restrict ourselves to the low-lying CEF states in the 
group s\. Within the T^-T2-T^ CEF model, the matrix ele- 
ments of J z = — 4cr, H" = 105£ and D S 4 = -630/7 are ex ~ 
pressed as 
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Fig. 2. (Color online) The mean-field p-T phase diagram. The dotted lines 
indicate the pure transition temperatures without other OPs. The boundary 
between the AFH and the AFM is of the first order. 





i cos 9 




-1 cos 8 


i sin 9 





-ism 8 




n = 



000 
i oj 



(4) 



(5) 



To be specific, we set the CEF energies as = 0, E% - 50 K 
and E(' = 170 K, which are roughly the same values as was 
used in the susceptibility analysis. 19 ' The mixing parameter 
of the CEF wave functions is given by 9 — 0.998. 

First, we discuss the AFH and the AFM mean-field solu- 
tions based on the nearest-neighbor exchange model under 
external magnetic field h, = -Agjp^H (gj = 4/5 in the LS- 
coupling scheme), 

H = 7 X + f Yj &j - h * Z at + £ H < CEF ' (6) 

<u> <<J> 
where z — 8 is the coordination number and Hf EF is the diag- 
onal matrix with (-Ej 1 ', E%, E^). The comment on the dotria- 
contapole r\ will be given later. 

Assuming the linear p dependences in J and D so as to re- 
produce the observed values of Th at ambient pressure, tri- 
critical point at p c = 1.5 GPa and the zero-T AFH- AFM 
boundary at p = 0.6 GPa, 11 ' we obtain the mean-field p- 
T phase diagram as shown in Fig. 2. Here, we have used 
J(p) = J c + a(p-p c ) and D(p) = D c +/3(p-p c )with J c = 40.6 
K, D c = 91.4 K, a = 1.27 K/GPa and /3 = 1.0 K/GPa. 

Since Fj 2) state is high enough as compared with the tran- 
sition temperatures, the multipoles in the P* 1 '^ states are re- 
duced to as <x — > cr x sin 9, £ — > cr y cos 9 (cr is the Pauli matrix) 
and i] is vanishing. H CE¥ — > cr z (E\ - Ei)j2 acts as the "mag- 
netic field" along z direction in the two low-lying CEF space. 
This situation is discussed in the Ginzburg-Landau framework 
by Haule and Kotliar. 26 ' Note that within the two low-lying 
CEF space the effective couplings at the tri-critical point have 
almost the same magnitude, i.e., J c sin 2 9 D c cos 2 9 with the 
pseudo U(\) (XY) symmetry. 



(3) 
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Fig. 1 . (Color online) The local wave functions under (a) the AFM order, (b) the AFH order, (c) the AF dotriacontapole order, and (d) the paramagnetic 
state. The shape and the color-map represent the charge and the magnetic-charge densities, respectively. 



In the AFH and the AFM phases, the local wave functions 
are represented by 



AFH) ~ llf^ifelTa), 
|AFM>~ \lf y )±a\T 2 ), 



(7) 
(8) 



respectively, where a is a real constant depending on the tem- 
perature and the model parameters and the upper and the 
lower signs correspond to A and B sublattices. Since |±4) oc 
exp (±i4<p) {if is the azimuthal angle) and II^) is pure imagi- 
nary, AFH) is a real function that breaks the mirror symmetry 
with respect to the plane parallel to the c-axis in the electric 
potential, while AFM) is a complex function representing a 
circular current with keeping the local symmetry in the elec- 
tric potential. In the j- j coupling scheme within j = 5/2 mul- 
tiplet, 48 ' these states are expressed by 

cosf?/ \ 
|AFH)~— \\d 5 )+3\d 3 ) + 2\di)] 



sin 



(l+4) + 1-4)) ± »3=(l+4> - 1-4)), 



(9) 



COS0/ \ 

|AFM) ~ -=(Ws> + 3 \d 3 ) + 2 \di)j 



sin 



Vl4 

(l+4) + 1-4)) ± ^(l+4> - 1-4)), (10) 

where |+4) = |+5/2) ® |+3/2), |-4) = |-3/2) ® |-5/2) and 
\d m ) = l+m/2) ® |— m/2). From these expressions, two or- 
dered phases are essentially stabilized by the overlap among 
|+4) + 1— 4) and |+4) - 1-4), and the relative phase in the linear 
combination discriminates two distinct orders. 

In Fig. 2, the dotted lines indicate the pure transition tem- 
peratures without other OPs, which indicates that those orders 
are meta-stable states slightly higher than the real ordered 
states. The lst-order thermodynamic boundary between the 
AFH and the AFM phases is determined by the comparison 
of the free energies. Note that the lst-order phase boundary 
has a strong sample dependence 8 ' 1011 ' maybe due to uncon- 
trollable inhomogeneity of samples under pressure. This is 
consistent with the fact that even in the AFH phase the AFM 




Fig. 3. (Color online) The mean-field p-H (H || c) phase diagram at T = 0. 



exists, whose volume fraction increases by applied pressure. 6 ' 
The suppression of the HO state due to the Rh impurities 50 ' 
is discussed in the view of the competing orders of the AFM 
and the AFH. 51 ' With use of the above J(p), D(p) and the 
CEF parameters, the p-H (H \\ c) phase diagram at T — is 
obtained as shown in Fig. 3. The p dependences of H H and 
H N are qualitatively agreement with the observed ones. 10 ' 

3.2 The inelastic neutron scattering spectra 

Next, we discuss the inelastic neutron scattering spectra 
of the collective excitations from the above mean-field solu- 
tions. A striking feature of the observed magnetic response is 
that both the CEF and the itinerant aspects are manifested. 4 ' 5 ' 
Namely, the sharp CEF-like excitations propagate along the 
tetragonal basal plane at low energies, whereas fluctuations 
propagating along the c-axis constitute the high-energy broad 
magnetic excitations. The dispersive excitations have a gap 
of 2 meV at the AF zone center and are strongly damped 
out above T a . In the HO phase, the well-defined propagating 
magnon develops. However, upon entering the AFM (7n = 
1 .5 K) by decreasing T at p = 0.67 GPa, neither quasi-elastic 
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Fig. 4. (Color online) The inelastic neutron scattering spectra at ambient 
pressure (in the AFH phase) in the nearest-neighbor exchange model. The 
dynamical structure function for the hexadecapole scattering is also shown 
(the dashed lines). The gray lines indicate the dispersion relations for the 
collective excitations. 
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Fig. 5. (Color online) The inelastic neutron scattering spectra at p = 1.2 
GPa (in the AFM phase). The dipole intensity is diminished near the or- 
dering vector Q = (1, 0, 0). 



nor inelastic responses at Q have been observed. 52 ' Thus, 
these excitations are characteristic of the HO phase. 10 ' 50 > 52 ' 

Using the generalized Holstein-Primakoff framework, 53-56 ' 
we calculate the dynamical structure function at T — 0, 

S(k,e) = —lmxz(k,e), (11) 
n 

where Xz(k> e ) ls trie dipole-dipole green function defined as 

Xz (k, e) = ((a; cr)) k (e) = -i dt ([oik, t), cr(-k)]) e ia . 

Jo 

(12) 

The dynamical structure function of the hexadecapole scat- 
tering is also calculated in a similar manner with g)) k (e), 
although it cannot be measured by present experimental tech- 
niques. 

Figure 4 shows the excitation spectra at T — in the AFH 
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Fig. 6. (Color online) The p dependence of the excitation energy and the 
intensities at T = 0. 



phase at ambient pressure along (100) direction. The intensity 
due to the magnetic scattering slightly increases toward the 
ordering vector Q, while the hexadecapole scattering inten- 
sity is rapidly decreasing toward Q. Similarly, the excitation 
spectra in the AFM phase at p — 1.2 GPa is shown in Fig. 5. 
In contrast to the AFH phase, the magnetic scattering inten- 
sity is almost diminished toward Q, while the hexadecapole 
scattering remains finite. In general, the collective mode is ex- 
cited by transverse components (with off-diagonal matrix ele- 
ments) with respect to the OP component. In the AFH phase, 
the magnetic dipole J z plays a role of the transverse compo- 
nent, while it is the hexadecapole that plays this role in the 
AFM phase. The latter is invisible by the neutron scattering. 
Here, it is again important that the low-lying CEF states do not 
contain the active dipole operators J x and J y . This mechanism 
for the missing intensity in the AFM phase is also pointed out 
by Haule and Kotliar, 26 ' and a similar discussion was given 
in the context of the doublet CEF model. 20 ' Note that the 
missing intensity in the AFM phase should be recovered away 
from Q as shown in Fig. 5. This provides a natural explana- 
tion of the mode observed in neutron scattering experiments 
under pressure. 4 ' 5 ' 10 ' 50 ' 52 ' 

The p dependence of the energy of the collective excita- 
tion at Q and its intensities is shown in Fig. 6. The increas- 
ing magnetic scattering intensity with increase of p suddenly 
drops by an order of magnitude across the lst-order transition 
to the AFM phase, at which the excitation energy (jl>q{Q) takes 
a minimum value. On the contrary, the invisible hexadecapole 
scattering is largely intensified across the phase boundary. 

Similarly, the H (|| c) dependence of the intensity is shown 
in Fig. 7. The excitation energy has a non-monotonous H de- 
pendence in the AFH phase. It is vanishing toward the second- 
order Hu since the magnitude of the OP is vanishing as well. 
The excitation energy gap is roughly given by the energy dif- 
ference of two ordered phases, and fluctuations around them 
constitute the pseudo Goldstone mode. 26 ' 

At this point, we shall mention the possibility of the dotri- 
acontapole order 77. Since the position of F (2) is unclear due 
to the lack of direct experimental information on the CEF 
scheme, one might consider that the AF 77 order would take 
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Fig. 7. (Color online) The H || c dependence of the excitation energy and 
the intensities at p = 0.8 GPa and T = 0. 



Fig. 8. (Color online) The H \\ c dependence of the ordered hexadecapole 
moment £ (Q) and the induced dipole and dotriacontapole moments. 



Table V. The angle dependence of the RXS amplitude in the nonrotated 
(ct-ct 1 ) and the rotated (tr-n') polarization channels both for the El and 
E2 transitions in the presence of the AFH order. The 3rd and the 4th rows 
indicated by (On) represent the angle dependence due to the induced AF 
O22 order under the in-plane magnetic field. 
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angle dependence 
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cos 2t^ 
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sin 8 sin 2l/i 
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(022) 0--JT' 



2V2 
2V2 



( 1 - cos 2ifr) 



sin 8 sin 2t^ 



(21 

over the AFH order, provided that the Fj state is low enough 
in energy. However, since 77 is connecting between two Ti 's 
while cr is connecting between V\ and F2 (see eqs. (3) and 
(5)), the occurrence of the AFM order has little influence on 
the 77 operator and vice versa. Hence the drastic change of the 
magnetic scattering intensity across the phase boundary does 
not occur in the case of the AF 77 order. This is confirmed by 
a similar calculation in the above (not shown). By this reason, 
the AF dotriacontapole order is not a plausible candidate for 
the HO phase. 

3.3 The RXS amplitude and the field-induced moments 

Let us now consider direct consequences of the AFH order. 
One of the powerful tools to observe the multipole order is 
the RXS measurement. 28 ' Since the multipole order accompa- 
nies a specific shape in the charge and/or the magnetic-charge 
densities different from the original lattice symmetry, the az- 
imuthal angle dependence of the RXS intensity provides a di- 
rect evidence for the OP symmetry. According to a concise 
formalism for the RXS amplitude, 29 ' which assumes a spheri- 



cal symmetry in the intermediate resonant processes of a core 
hole, the angle dependences of the scattering amplitude are 
summarized in Table V. Here, 6 is the scattering angle and 
if/ is the azimuthal angle, whose origin is chosen when the re- 
polarization vector e a is parallel to a axis in the case of the ro- 
tation axis Q = (001). The detailed experimental geometry of 
the scattering is shown in Fig. 1 of ref. 29. Note that although 
Q = (001) and (100) are essentially equivalent momenta in 
the Brillouin zone, the azimuthal dependences with respect to 
these two rotation axes are different. The 3rd and 4th rows in 
Table V indicated by {O22) will be explained below. 

The presence of the AFH order causes an additional elec- 
tric field which breaks the D4/, symmetry at U-site. 34) The ad- 
ditional electric potential at U-site is given by 



(13) 



where (r 4 ) is the average of r 4 over the 5/ radial wave func- 
tion. In the LS -coupling scheme, c - -14 V35/363. The re- 
sultant electric potentials generated at Si- and Ru-sites are 
proportional to xy(x 2 -y 2 ) and xy, respectively. This is consis- 
tent with the space group, No. 128 {PA jmnc) corresponding to 
the AFH order. 32 ' The 4-fold axis at Ru-sites is lost by the xy- 
type electric potential, which should result in a finite asymme- 
try parameter (77) in the electric field gradient tensor. However, 
the NQR frequency at Ru-site has shown no changes at the 
onset of the HO within experimental accuracy. 33 ' The lack of 
the frequency shift may be ascribed to smallness of the elec- 
tric field due to the higher-rank hexadecapole order. In fact, 
the ratio of <pti{r) and the electric potential from the 2nd-rank 
(O20) in the CEF potential, fair), is roughly estimated as 



(>-) 



Mr) 



10~ 3 x 



r 2(yiy 



(14) 



which is at least three order of magnitude small. This small- 
ness is also consistent with the fact that no measurable lattice 
distortions have been observed. 

When the magnetic field along the c-axis is applied, the AF 
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Fig. 9. (Color online) The in-plane magnetic-field angle dependence of the 
uniform ordered moments. 
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Fig. 10. (Color online) The in-plane magnetic-field angle dependence of 
the staggered ordered moments. 



dotriacontapole is induced through the 3rd-order term, 



(15) 



in the Ginzburg-Landau (GL) free energy. The H (|| c) depen- 
dence of cr(0), £(0) and 77(g) at T = are shown in Fig. 8. 
The induced AF dotriacontapole has a maximum around the 
field at which cr(0) x £(0) becomes the largest. The induced 
rj(Q) causes an additional internal magnetic field along c-axis 
at Si-site, but not at Ru-site. It is interesting to note that the 
single-crystal 29 Si NMR line shows an additional broadening 
in the temperature range [14, 17.5] K, 57) which may be rele- 
vant to the induced rj(Q). 

We consider the effect of a magnetic field in the basal plane. 
Since the GL free energy contains the 4th-order coupling, 



(16) 



the in-plane magnetic field should induce the AF O22 (x 2 -y 2 )- 
type quadrupole order. However, the induced moments are 
expected to be small since they are almost inactive at low 
temperatures. In order to elucidate the induced orders semi- 



quantitatively, we adopt the mean-field theory with the full 
CEF states under eq. (1). The in-plain field angle dependences 
[H = //(cos if, sin if, 0)] of the uniform and the staggered 
components of J x , J y and O22 are shown in Figs. 9 and 10. 
Here, we have used H — 5 T at T = IK, and the CEF parame- 
ters, B 20 = -5.18, B40 = -0.744, B 44 = 0.546, B 60 = B M = 0, 
which reproduce the low-lying CEF states in the T^^-Tz-T?^ 
model. All of the induced dipole moments, J X (Q) and J V (Q), 
are extremely small because of their inactiveness at low tem- 
peratures. The induced AF O22 quadrupole moment becomes 
the largest for the applied magnetic field along [110]. Note 
that the uniform quadrupole 022(0) is also induced, although 
it is quite small. It should be stressed that the induced AF 
quadrupole 022(6) i s a l so detectable by the RXS. It is worth 
performing the RXS under the in-plane field along [110] even 
in case that the RXS intensity of the AFH itself is too small 
to detect. The angle dependence of the scattering amplitude is 
summarized in the rows indicated by (O22) in Table V. The in- 
plane magnetic-field measurements on the elastic response in 
the transverse (en -Cn)/2 mode and the thermal expansion (it 
has already been performed without magnetic fields, 58 ' show- 
ing the intrinsic softening in this mode) are highly desirable. 

Similarly, a uniaxial stress can also induce AF quadrupoles. 
Namely, through the 3rd-order coupling in the GL free energy, 

T, ~ 0^(0^(6)022(0), or r 3 ~ 0^(0)^(6)022(0), (17) 

the AF quadrupoles of the (x 2 - y 2 )-type and the (xy)-type 
are induced by applying the uniaxial stresses along [110] and 
[100] directions, which correspond to the strain field of e xy 
and e v = e xx - e yy , respectively. Thus, a detection of these in- 
duced AF quadrupole under the uniaxial stress also provides a 
useful test for the AF hexadecapole order, although the mag- 
nitude of the induced quadrupole moment may be very small 
due to inactiveness of these multipoles in low temperatures. 

The presence of the induced AF 022 quadrupole breaks 
the 4-fold symmetry at U-site. Quite recently, the emer- 
gence of two-fold oscillations in the magnetic torque un- 
der in-plane field rotation are sensitively detected in small 
pure crystals. 59 60 ' It should be stressed that the induced 022 
quadrupole with the two-fold symmetry does not necessar- 
ily result in the two-fold oscillation in the magnetic torque, 
since the induced magnetic and quadrupole moments them- 
selves depend on the magnetic field direction. Nevertheless, 
a relevance between the two-fold oscillations and the induced 
022 quadrupole might be interesting to address in future in- 
vestigations. 

4. Summary and Outlook 

We have discussed the nature of the AF hexadecapole order 
and its consequences on the basis of the localized /-electron 
exchange model. The hexadecapole moment has been se- 
lected as a plausible candidate for the hidden order parameter 
in view from the consistency with the strong anisotropy in the 
magnetic susceptibility and the magnetization process. This 
consistency requires the large energy separation between the 
low-lying CEF singlets and the excited doublets, which ex- 
cludes all the multipoles in two-dimensional representations 
from candidates for the hidden order parameter. 

The mean-field and the collective-excitation analyses pro- 
vide reasonable agreements with the characteristic of the 
phase diagram and the missing neutron-scattering intensity 
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out of the hidden order phase. No significant change in the 
Ru-NQR frequency upon entering the hidden order phase can 
be understood by smallness of the electric potential associ- 
ated with this higher-rank order parameter. Nevertheless, this 
statement is unambiguously justified by more elaborate the- 
oretical and/or experimental arguments on the electric poten- 
tials quantitatively. The magnetic fields along the c-axis in- 
duces the AF dotriacontapole moment, which brings about 
additional internal magnetic fields at Si-sites of the same sym- 
metry as the dipole ones. Similarly, the AF O22 quadrupole 
moment should be induced by the in-plane magnetic fields, 
which affects the thermal expansion and the elastic constant 
of the transverse (c\ \ -c\i)j2 mode. The uniaxial stress along 
[110] and [100] directions can also induce the AF On and 
O xy quadrupoles, respectively. It would be interesting to de- 
tect these induced AF quadrupoles by the RXS under the in- 
plane magnetic field or the uniaxial stress. 

From general group-theoretical point of view, it is allowed 
a possibility that an order parameter in the high- p AFM phase 
belongs to the subgroup of No. 128 space group correspond- 
ing to the low-p AFH phase. In this case, the high-/? phase is 
characterized by the coexistence of the AFM and the AFH, 
and the phase transition between the low-/? and the high-/? 
phases can be the second order. Although the present work 
based on the localized exchange model does not have such a 
coexistent-phase solution, this possibility can provide a natu- 
ral explanation for the similarity of the Fermi surfaces across 
the pressure-induced phase transition. 

Since the superconductivity does not coexist with the 
AFM phase, it is presumable that the pairing is mediated by 
the pseudo Goldstone mode via the dipole coupling in the 
non-magnetic hexadecapole ordered state. 61-63 ' The marked 
anisotropy in the upper critical field H C 2 in spite of rather 
3 -dimensional character of fermi surfaces 1 ' 64 ' can reconcile 
with the scenario of the magnetic exciton mediated supercon- 
ductivity, which is sensitively damped against the magnetic 
field along c axis but is insensitive for the in-plane fields (see 
Fig. 7 for example). Thus, it would be interesting to elucidate 
an interplay between the magnetic excitations and the super- 
conductivity in future investigations. 
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